Possible Formation Scenario of the Quark Star of Maximum Mass 

around 0.7 M 
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ABSTRACT 

If there exists the quark star of maximum mass ~ 0.7Mq as suggested by recent Chandra 
observations, we show that the general relativistic collapse of a neutron star of mass greater than 
the maximum mass of the neutron star with angular momentum parameter q = cJ/GM > 1 
may lead to such a strange star. Here J and M are the angular momentum and the gravitational 
mass of the neutron star, respectively. Under the cosmic censorship hypothesis, such a star can 
not be a black hole directly. The jet formed in the soft core might explode the outer envelope 
and leave the quark star of mass ~ 0.7M Q . The remnant quark star has < 10 53 erg rotational 
energy so that the formation of the quark star may be related to the central engine of GRBs. 
The detailed numerical simulations are urgent to confirm or refute this scenario. ( version 3 May 
2002) 

Subject headings: stars: neutron 
1. Introduction 

Recently Drake et al. (2002) have made Deep 
Chandra LETG+HRC-S observations of the iso- 
lated neutron star candidate RXJ 185635-3754 and 
found that the X-ray spectrum is well-represented 
by a ~ 60eV black body. Their data contain 
no evidence for pulsation. The derived column 
density favors the distance to RXJ185635-3754 
to be D ~ 140pc instead of D ~ 60pc (Walter 
2001). Although Walter & Lattimer (2002) claim 
that the distance is D ~ 120pc and the heavy- 
element atmosphere of Pons et al. (2002) should 
be taken, we here follow the arguments by Drake 
et al. (2002). 

From the observed flux, the distance and 
the observed temperature of RXJ185635-3754, 
only the radiation radius defined by = 
R/y/l — Rg/R is determined since the observed 
luminosity and the observed temperature suffer 
from the gravitational redshift. For the given i?oo, 
M is given as a function of R as 



maximum value R^/GSy/S = 0.7 M Q (R 00 /5Mm) 
at R = Roo/\/S. If Roc is ~ 5.6kmZ?i 4 o as sug- 
gested by Drake et al. (2002), the maximum 
mass is ~ 0.7MqDhq where D140 is D in the unit 
of 140pc. Therefore Roo can not determine the 
radius of the star but gives us the upper limit of 
the mass as well as the upper limit of the radius 
of RXJ185635-3754. 

X-ray luminosity of RXJ185635-3754 is - 6 x 
10 31 erg/sD 2 40 . If this is due to the accretion from 
a uniform medium (Bondi & Hoyle 1944) we have 
the relation as 
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where V is the velocity of RXJ185635-3754 and e 
is the conversion efficiency of the accreting matter 
to the radiation and is ~ 0.3c 2 for the quark star 
with the maximum mass ~ 0.7M©(Witten 1984). 
From the above equation we have 
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M is positive only for < R < R^ and has the 
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Since the proper motion of RXJ 185635-3754 is 
0.3 arc-second/y - 200Z)i 40 kms- 1 (Walter 2001) 
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the required number density of the interstellar 
matter is rather high even for M ~ Mq . While 
the column density to RXJ 185635-3754 is esti- 
mated as ~ 10 20 cm- 2 (Drake et al. 2002) so 
that the high density region should be local with 
the maximum size of ~ 100AU. Such a local high 
density region of the size < 100AU is observed 
in the interstellar medium using 21cm HI line by 
Frail et al. (1994) so that RXJ185635-3754 may 
be passing through the high density region now. 

If M is <~ 0.1M Q , the required density of the 
interstellar matter is so high that the column den- 
sity may be in conflict with the observed one or we 
need very small ( ~AU) high density region in the 
interstellar matter. In the latter case we expect 
the change of the luminosity in a time scale of 10 
days or so, while only ~2% change of the luminos- 
ity has been observed in the 19 months (Drake et 
al. 2002) so that the mass is much greater than 
<~ 0.1M Q . In this Letter, assuming the mass and 
the radius of RXJ185635-3754 are ~ O.7M and 
<~ 3.2km, we propose a possible formation scenario 
of such an object. 

2. Spherically Symmetric Model 

First such a small radius object can not be the 
neutron star ( Lattimer & Prakash 2000). Even for 
the softest equation of state of the high density nu- 
clear matter, R^, is larger than 10km ( Lattimer & 
Prakash 2000). Except for the free neutron case ( 
Oppcnhcimcr & Volkoff 1939), the maximum mass 
of the neutron star is greater than 0.7M Q ( Lat- 
timer & Prakash 2000), which is in conflict with 
the maximum mass of RXJ185635-3754 from the 
observed i?^. The strange quark star is one of the 
candidate for such an object like RXJ185635-3754. 
In the quark star the baryon number is distributed 
among three quarks u, d and s so that the degener- 
ate energy decreases by a factor of 0.89 compared 
with the proton and neutron nuclear matter (Wit- 
ten 1984). In this sense the strange quark mat- 
ter is stable although it contains unstable strange 
quarks. This reason is essentially the same as that 
the unstable neutron exists stably in the stable nu- 

1 Here we assume that the optical counterpart is the same 
object as RXJ185635-3754 although the luminosity in the 
optical band can not be explained by the Rayleigh-Jeans 
tail of ~ 60eV black body spectrum of RXJ185635-3754 
(Pons et al. 2002; Drake et al. 2002; Walter & Lattimer 
2002). 



cleus. In MIT bag model with the large bag con- 
stant B = (245MeV) 4 suggested from the finite- 
temperature QCD lattice data (Peshier, Kampfcr 
& Soff 2001) , the maximum mass of the quark 
star Mmax is given by 

M - = °- 7M ^(2l5W rl/2 ' (4) 

R = :! - !,/ '-''' : -2T^' " < B > 

Pc = l-5xl0V m - 3 (p^),(6) 

where R and p c arc the radius and the central 
density of the maximum mass quark star (Wit- 
ten 1984). Note here that the central density 
of the quark star with Mmax ~ 0.7M Q is much 
higher than the central density of the neutron star 
<~ 10 15 gcm~ 3 . Therefore to make a quark star 
from the neutron star we must consider how to 
increase the central density of the neutron star. 
If the bag constant is smaller, for example B = 
(145MeV) 4 as usual , the maximum mass and the 
radius are 2M Q and 11km, respectively. However, 
in this case with R x =5. 6km , M becomes small 
~ 0.1 5M Q which is in conflict with the X-ray lumi- 
nosity of RXJ185635-3754 as stated in §2 so that 
we prefer the larger bag constant. 

A naive way to increase the central density of 
the neutron star is to add the mass to the max- 
imum mass neutron star by accretion from the 
companion star or the fall back gas in the su- 
pernova event ( MacFadyen , Woosley , & Hegcr 
2001). However the maximum mass of the neu- 
tron star is 1.4M ~ 2M ( Lattimer & Prakash 
2000), which is at least two times larger than the 
putative maximum mass of the quark star rele- 
vant to RXJ185635-3754. The most plausible out- 
come of such a collapse is not the formation of 
the quark star but the black hole. In a sense the 
soft equation of state of the strange quark matter 
(P = (p — 4£?)/3) is favorable to the formation 
of the strong shock in the supernova explosion as 
shown by Takahara & Sato (1988). If the equation 
of state is soft, the large density gradient is accu- 
mulated to produce the strong shock while for the 
hard equation of state, the small density gradient 
is enough to produce the pressure gradient to sup- 
port the falling gas in the collapse. However the 
shock is the strongest just before the formation of 
the black hole, that is, when the mass of the col- 
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lapsing star is just a bit smaller than the maximum 
mass of the neutron star in the spherically sym- 
metric model (Van Riper & Arnett 1978). Since 
in our case of the formation of the quark star, the 
mass of the collapsing matter is much larger than 
the maximum mass of the quark star so that wc 
expect the formation of the black hole. The de- 
tailed numerical simulations might be needed to 
confirm that the shock can not explode the outer 
envelope and the black hole is formed. 

3. Rotating Model 

I propose another promising scenario. First 
let us review the simulations of general relativis- 
tic collapse of rotating stars. Nakamura (1981) 
first performed the numerical simulations of gen- 
eral relativistic collapse of axially symmetric stars 
of mass 10M Q leading to the formation of rotat- 
ing black holes. In these simulations the equa- 
tion of state is hard with 7 = 2 in the limit of 
P 3> Pnuc and the rotation is differential with var- 
ious values of initial angular momentum (Naka- 
mura 1981). Nakamura & Sato (1981a) performed 
similar simulations with soft equation of state of 
7 = 4/3 for two kinds of rotation laws with various 
values of angular momentum. The results of these 
simulations crucially depend on one parameter q 
defined by 

where J and M are the total angular momen- 
tum and the gravitational mass of the initial ro- 
tating star, respectively. When q < 1 the black 
hole is formed irrespective of initial rotation laws, 
the equation of state and the initial density dis- 
tribution (Nakamura & Sato 1981b). If q > 1 , 
the black hole is not formed but the expanding 
ring, the expanding disk or the jets appear de- 
pending on initial rotation laws, equation of state 
and initial density distribution. This conclusion 
was confirmed by Stark & Piran (1985) and Shi- 
bata (2000). 

Physical reason for the above result is as fol- 
lows; Let us consider the collapsing star of mass 
M and the angular momentum J. The radius Rb 
when the centrifugal force balances the gravita- 
tional force is given by 



Using q, we can express Rb as 
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Equation (9) means that if q < 1, Rb is smaller 
than the black hole, which suggests that when the 
rotation becomes important, the star itself is al- 
ready inside the black hole. Therefore it is natural 
that the black hole is formed irrespective of initial 
rotation laws, the equation of state and the initial 
density distribution. 

In the sense of general relativity, the results of 
the simulations support the cosmic censorship hy- 
pothesis (Penrose 1973) which says that the sin- 
gularity does not exist outside the event horizon 
in nature. Under this hypothesis, the Kerr metric 
with a/m = q < 1 is known to be unique as the 
ultimate space-time of the gravitational collapse 
(Robinson 1975). Note here that q defined in Eq. 
(7) and a/m in Kerr metric are the same quan- 
tity. If the whole of the collapsing star with q > 1 
becomes a black hole, that is the counter exam- 
ple of the cosmic censorship hypothesis. However 
numerical simulations have not yet found such a 
case. 

In relation to the formation of the quark star, 
Nakamura & Sato (1981a,b) showed an interesting 
case of the jet formation. In their model of A146 
which is the general relativistic collapse of the ro- 
tating star with 7 = 4/3 equation of state and the 
initial q = 1.46, the jet appeared in the final phase 
of the collapse. This is in accord with the result 
of Van Riper & Arnett (1978); Takahara & Sato 
(1988). In the rotating collapse, when the size of 
the star becomes ~ Rb in Eq. (9), the collapse 
perpendicular to the rotation axis halts while the 
collapse along the rotation axis proceeds. When 
the equation of state is soft, the large density gra- 
dient accumulates to produce the jet. While if the 
equation of state is hard, this is not so. In reality 
Nakamura, Oohara & Kojima (1987) performed 
the simulations of the accreting neutron star with 
rotation. They used the equation of state with the 
maximum mass 1.4M Q and simulated the collapse 
of the rotating neutron star of mass 1.9M Q . They 
found no jets. The final result was the differen- 
tially rotating neutron star. 

Recently Baumgarte, Shapiro & Shibata (2000) 
computed the equilibria of differentially rotating 
neutron stars. They stated that some of the mod- 
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els exceed the Kerr limit, that is , q > 1. If we add 
certain amount of matter to exceed the maximum 
mass neutron star keeping q > 1, the collapse will 
start. At Rb, the density is estimated as 



(^nugget) should be 



Pb 



M 



- 6.7 x 10 16 gcm- 3 / -1 ^ 



-6/ M x-2 



(10) 

where / <~ 1 is the form factor of the rotating 
star. Since p^, can be in the quark phase for ap- 
propriate parameters, the equation of state in the 
central core can be soft at the centrifugal bounce. 
Then the matter will collapse along the rotational 
axis to form the jet like A 146 in Nakamura & 
Sato (1981a, b). If the jet is strong enough to 
explode the outer envelope and leave the central 
core of mass ~ O.7M , the strange quark star 
might be formed. The important point here is 
that if the cosmic censorship hypothesis is correct 
as suggested by numerical simulations, the whole 
of the rotating neutron star of mass greater than 
the maximum mass of the neutron star with q > 1 
can not be a black hole directly while if q = it is 
possible to be a black hole directly. 

The detailed numerical simulations are urgent 
to confirm or refute the present scenario of forma- 
tion of the strange quark star of maximum mass 
around 0.7 M . 

4. Discussion 

Olinto (1987) showed that a seed of strange 
matter in a neutron star will convert the star into 
a strange star. As a seed, the quark nugget formed 
at the QCD phase transition in the early universe 
(Witten 1984) can be considered. Then the cap- 
ture cross section of the quark nugget by the neu- 
tron star is given by 
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where i?NS> A^ns and V re i are the radius of the 
neutrons star, the mass of the neutron star and 
the relative velocity, respectively. In order that 
at least a single quark nugget should be captured 
in the age of the universe {tjj) for a certain neu- 
tron star , the number density of the quark nugget 
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Since the mass density of the quark nugget is 
at most the same as the dark matter density, 
each mass of the quark nugget m nU ggot should be 
smaller than ~ 10 12 g. However such a small mass 
quark nugget may not be formed due to the evap- 
oration (m nugg ot > 10 20 g for survival; Bhattachar- 
jee et al. (1993)). Therefore to explain the pos- 
sible strange quark star of RXJ185635-3754, the 
quark nugget should have been captured by chance 
by the progenitor neutron star. 

The putative mass of <~ 0.7M Q is just the 
mass of the typical white dwarf (Bragaglia Renzi 
& Bergeron 1995). Therefore if wc can com- 
press the white dwarf of density ~ 10 7 gcm -3 to 
~ 10 16 gcm~ 3 , the strange quark star of mass 
O.7M might be formed. A possible process is 
the tidal pinching proposed by Luminet & Pichon 
(1989). If the moderate massive black hole of mass 
~ 2000M Q exists, a white dwarf passing inside the 
tidal radius of the black hole may be pinched and 
the density increases by more than a factor 50 ( 
Luminet & Pichon 1989). Although the existence 
of such an intermediate mass black hole has been 
suggested (Matsumoto et al. 2002), the number 
density of the white dwarfs is small ~ O.Olpc -3 so 
that the pinching number in the age of the uni- 
verse for a single intermediate mass black hole is 
very small as 

QGMbh GMbh . 

pinch ~ n>WD( j A — T72 ) V rcltu 
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Therefore to explain the possible strange quark 
star of RXJ185635-3754, the white dwarf should 
have been tidally pinched by the intermediate 
mass black hole by chance. 

In any scenario, if the quark star of mass M qs <~ 
O.7M is formed, the huge binding energy ~4x 
10 53 erg is liberated. The major part of this energy 
may be used to explode the outer envelope of the 
progenitor neutron star. If , however, only ~ 10% 
of this energy is deposited as the kinetic energy, 
it is enough to explain the luminous supernova 
such as SN1998bw. In the rotating model, the 
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remnant quark star should have an initial value of 
q qs < 1. Then the rotational energy is estimated 
as E rot ~ {GM qs /Rc 2 fq 2 qs M qs c 2 ~ 10 53 erg. This 
energy is larger than or comparable to the energy 
needed for GRB so that the formation of the quark 
star may be related to the central engine of GRB. 
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